Fast, Runaway Evaporative Cooling to Bose-Einstein Condensation in Optical Traps 
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We demonstrate a simple scheme to achieve fast, runaway evaporative cooling of optically trapped 
atoms by tilting the optical potential with a magnetic field gradient. Runaway evaporation is possible 
in this trap geometry due to the weak dependence of vibration frequencies on trap depth, which 
preserves atomic density during the evaporation process. Using this scheme, we show that Bose- 
Einstein condensation with ~ 10 5 cesium atoms can be realized in 2 ^ 4 s of forced evaporation. 
The evaporation speed and energetics are consistent with the three-dimensional evaporation picture, 
despite the fact that atoms can only leave the trap in the direction of tilt. 



PACS numbers: 67.85.Hj, 64.70.fm, 67.85.-d 

The possibility to manipulate Bose-Einstein conden- 
sates (BECs) and degenerate Fermi gases of cold atoms 
in optical traps opens up a wide variety of exciting re- 
search; prominent examples include spinor condensates 
[l| , Feshbach resonance in cold collisions Q , and BECs of 
molecules 0, H|. In many early experiments, condensates 
were first created in a magnetic trap and subsequently 
transferred to an optical dipole trap. These experiments 
could be greatly simplified after direct evaporation to 
BEC in optical traps was demonstrated In this paper, 
we describe a further improvement on dipole-trap based 
evaporation, which allows for runaway cooling without 
significant increase in trap complexity. 

Evaporative cooling proceeds by lowering the depth 
of a confining potential, which allows atoms with high 
kinetic energy to escape and the remaining particles to 
acquire a lower temperature and higher phase space den- 
sity through rethermalization. Starting from a sample 
of precooled atoms in a dipole trap, one can in principle 
perform forced evaporative cooling on optically trapped 
atoms by constantly reducing the trap depth until quan- 
tum degeneracy is reached. This method has been suc- 
cessful in creating rubidium BEC in a dipole trap, and 
has become a critical component in recent experiments 
on quantum gases of Cs Q, Li Q, K Q and Yb @. 
In all these experiments, forced evaporative cooling in 
the dipole trap is realized by reducing the intensity of 
the trapping beam, and consequently also the restoring 
forces. In later discussion, we will refer to this approach 
as trap- weakening scheme. 

Evaporative cooling in optical traps remains one of the 
most time-consuming and technically challenging steps in 
condensate production. Fundamentally, this is due to the 
fact that cooling by weakening the trapping potential in- 
evitably reduces the collision rate. Here runaway (accel- 
erating) evaporation is essentially impossible even with 
perfect evaporation efficiency and purely elastic collisions 
[lQj . Within experimentally accessible times, the trap- 
weakening method puts a severe limit on the maximum 
gain in phase space density one can reach. Several auxil- 
iary schemes have been successfully implemented in order 
to increase the evaporation speed, including the dimple 
trap [6] and a zoom lens system [Tl| . These methods of- 
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FIG. 1: (color online) Trap-tilt based evaporation and ex- 
perimental apparatus, (a) Trap depth U decreases when an 
external potential gradient is applied to the optically trapped 
atoms, (b) Apparatus for evaporation of cesium atoms (black 
dot) in a crossed-beam dipole trap. A strong, slowly- varying 
magnetic field gradient B'(t) over- levitates the atoms with 
magnetic moment \i against gravitational pull mg and evap- 
orates them upward. 



ten increase the complexity of the apparatus or require 
delicate optical alignment or manipulation. 

In this paper, we report a new and simple evaporative 
cooling scheme which can be immediately implemented 
in many existing experiments. Instead of reducing the 
intensity of the trapping beam, we reduce the trap depth 
by applying an external force on the optically trapped 
atoms. This trap-tilting method entails only a weak re- 
duction in confinement strength over a large range of po- 
tential depth and can significantly speed up the cooling 
process. Using this method, we demonstrate runaway 
evaporative cooling in a large volume dipole trap and 
reach Bose-Einstein condensation of cesium significantly 
faster than previous results [12[ . Finally, we comment on 
the conditions for runaway evaporation in a tilted trap 
and investigate the dimensionality of atomic energy se- 
lection in the evaporation. 

For this study, cesium atoms are first slowed by a Zee- 
man slower, collected in a magneto-optical trap (MOT) 
for 2 s, molasses precooled, and finally cooled and spin 
polarized by degenerate Raman-sideband cooling (dRSC) 
|13( to the lowest hyperfine ground state \F = 3,rriF = 
3), where F is the total angular momentum and is 
the magnetic quantum number; the apparatus for dRSC 
follows that in [13[. A crossed dipole trap and magnetic 
field gradient are employed to levitate and collect the 
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FIG. 2: Performance of trap-tilting based forced evaporation: 
(a) phase space density, (b) collision rate, (c) particle number 
and (d) density profile (d). Two evaporation paths: 4 s (solid 
dots) and 1.8 s (open circles) are shown. The dashed line in 
(a) shows simple exponential increase. In (d), time-of- flight 
absorbtion images and single-line optical density profiles are 
taken from the 1.8 s evaporation path. The expansion time 
is 70 ms, and the field of view is 1.2 mm x 1.2 mm. 



cooled atoms. The dipole trap is formed by intersect- 
ing two laser beams on the horizontal (x — y) plane; both 
beams are extracted from a single- mode, single frequency 
Yb fiber laser operating at the wavelength of 1064 nm, 
frequency offset by 80 MHz, focused to a 1/e 2 beam diam- 
eter of 540 fim (620 /im) and intensity of 1.9 W (1.6 W) 
in the y—(x—) direction. In the absence of trap tilt, 
the trapping frequencies near the bottom of the poten- 
tial well are = 2tt x (17,34,38) Hz. Dur- 
ing the dipole trap loading process, we switch on a uni- 
form magnetic field of 58 G in the (vertical) ^-direction 
to improve the atom number following the loading pro- 
cess [6] and apply a levitating magnetic field gradient 
of B' c = mg/ [1=31.3 G/cm, where mg is the gravita- 
tional force, \i = 0.75 \±b is the magnetic moment of 
the atoms in |3,3), and hb is Bohr magneton. After 
1 s of thermalization and self-evaporation in the dipole 
trap, we ramp the magnetic field to 20.8 G, where three- 
body loss is minimized [l4|, and the field gradient to 
37.8 G/cm in 1.85 s and begin our study on forced evap- 
oration. At this point, which we define as time t = 0, 
there are No = 1.9 x 10 6 atoms in the trap with a 
temperature of To = 470 nK, peak atomic density of 
no = 3.8 x 10 12 cm -3 , peak collision rate of Tq =133 /s. 



The background collision rate are below 1/60 s. 

We perform forced evaporative cooling by linearly in- 
creasing the magnetic field gradient B' from 37.8 to 
41.4 G/cm in 2.2 s and then to 43.5 G/cm in another 
3 s, which reduces the trap depth from 3.0 jaK to 1.0 /iK 
and then to 170 nK, as calculated from our potential 
model. The magnetic field and dipole trap intensity are 
kept constant throughout the process. To evaluate the 
cooling performance, we interrupt the evaporation at var- 
ious times to measure the particle number N, temper- 
ature T and trap frequencies w x ,y,z- Particle number 
and temperature are extracted from absorption images 
taken at low magnetic fields, following a 70 ms time- 
of- flight expansion at B = 17 G to minimize the col- 
lisions and B'{z) = B' c to levitate the atoms. Trap 
frequencies are measured from small amplitude oscilla- 
tions of the atomic momentum by abruptly displacing the 
trap center. Peak phase space density is calculated from 
(j) = n\^ dBl where n = Nco x co y u; z (mXdB/h) 3 is the peak 
atomic density, Xjb — h(27rmkBT)~ 1 ^ 2 is the thermal de 
Broglie wavelength, fc# is the Boltzmann constant and 
h is the Planck constant. Collision rates are calculated 
as T = n(av), where the elastic collision cross section is 
a = 87ra 2 , scattering length at 20.8 G is a = 200 ao [Hj], 
(v) = (WkBT/nm) 1 / 2 is the mean relative velocity. 

After 4 s forced evaporative cooling, we observe Bose- 
Einstein condensation from the appearance of bimodality 
and anisotropic expansion in time-of- flight images. At 
this point, the temperature is 64 nK and total particle 
number is 5 x 10 5 . An almost pure condensate with 10 5 
atoms was obtained after another 2.5 s. In this evapora- 
tion process, the mean truncation parameter is calculated 
to be fj = (JJ/ksT) = 6.5(3), the evaporation efficiency 
is Jev = — log(0/0o)/ log(A^/A^o) = 3.4. We observe an 
increasing collision rate and accelerating evaporation, in- 
dicating achievement of runaway evaporation; see Fig. 2. 

An alternative evaporation path is developed to mini- 
mize the time to reach Cs BEC. After a shorter magnetic 
field ramping process of 1 s, we ramp the field gradient 
from 38.9 G/cm at t = to 41.3 G/cm in 0.5 s and then 
to 43.5 G/cm in another 1.5 s. Here we reach BEC in 
as short a period as 1.8 s forced evaporation. Another 
1 s evaporation allows us to obtain 4xl0 4 atoms in an 
almost pure condensate, see Fig. 2 (d). Despite the rapid 
increase of phase space density, the collision rate actually 
decreases by 25% at the end of evaporation. The trun- 
cation parameter and evaporation efficiency are fj = 4.6 
and j ev = 1.9, respectively. 

Throughout both evaporation processes, the peak den- 
sity is moderate, n < 1.5 x 10 13 cm -3 . The collision loss 
rate, determined from the three-body recombination pro- 
cess [HI, is below 1/40 s at 20.8 G. Trap loss from colli- 
sions is negligible in the following discussion. 

To understand the advantage of the trap-tilting 
scheme, we analyze evaporative cooling in a model poten- 
tial. We combine the magnetic gradient potential and the 
gravitational potential as — jmgz, where 7 = B' / ' B' c — 1. 
The total potential V(x, y, z) can be modeled as 
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FIG. 3: Depth and oscillation frequency of a tilted trap, (a) 
shows the calculated normalized trap depth and frequencies 
(jj z and uj x = ujy as a function of the tilt £, based on Eq. (1). 
In (b), mean trap frequencies are plotted against the trap 
depth for a tilted trap (solid dots) and for a weakened trap 
(dotted line). The solid line shows a polynomial fit to the 
mean frequency, see Eq. (2). 




Evaporation time (s) 

FIG. 4: Evaporation speed: experiment (4 s path, circle dot) 
and models. We assume an initial collision rate of ro=133 /s, 
77 = 77 = 6. 2^6. 8 and no collision loss. Shaded area covers 
ID evaporation region with < v < 1 and all possible 77. 



V = -^[ e - 2 (- 2 +* 2 )A" 2 + e -2(y 2 +* 2 )/- 2 ] _ imgZj (1) 

where the first two terms come from the two horizontal 
trapping beams, the last term is the tilt potential. Here, 
we assume the two beams have the same beam waist w 
and peak light shift Uq/2 for convenience. 

We introduce the tilt parameter ( = e 1 / 2 jmgw/2Uo 
to parameterize the trap depth U and trap frequencies 
w x ,y,z- Using Eq. (1), the trap depth and frequencies 
are evaluated as a function of as shown in Fig. 3 (a). 
All quantities are normalized to those of an untilted po- 
tential, where the trap depth is Uq : and the tr ap fre- 
quencies (j Q z = V^oj® = V2cUy = y/4Uo/mw 2 . Note 
that the trap is unstable when £ > 1. In the range of 
10 -3 < U/Uq < 1, the geometric mean of the trap fre- 
quencies Q = (uOxUUyUUz) 1 ^ varies with the trap depth 
approximately as, see Fig. 3 (b), 

cD/cJ «1.05([//f/o) - 075(1) , (2) 

where uJq is the mean frequency of an untilted trap. 

The key to fast, runaway evaporation in a tilted trap 
lies in the gentle, almost negligible weakening of the trap 
confinement when the trap depth decreases. As the trap 
depth reduces by a factor of 100, the trap frequency only 
decreases by 45% in the z— direction and 14% in the 
other two directions. This should be contrasted to the 
trap-weakening method, which reduces trap frequencies 
by a factor of 10 under the same condition. In general, a 
weakening trap with O cx U v and v = 0.5 shows a much 
stronger dependence on the trap depth than the tilting 
trap with v = 0.075. 

Collision rate in a harmonic trap depends on the parti- 
cle number, trap frequencies and temperature. Assuming 
77 = U/ksT > 6 is a constant, we have 

T cx Nuj^T' 1 

cx U^U^U- 1 = (3) 



where (3 = 1/a + 3u — 1, a > parameterizes the tem- 
perature decrease by evaporating one atom [l6[ , namely, 

dlogT _ 77 + ft- 3 
a ~dI^N~ 3-3i/ ' [ ) 

and n > depends on the dimension of evaporation as 
discussed below. 

The condition for runaway evaporation is given by 
(3 < 0. For the trap weakening scheme with v — 1/2, (3 is 
positive for all 77. Runaway evaporation is thus impossi- 
ble. For the tilting scheme with v = 0.075, the exponent 
/3 is negative when 77 + k > 6.58, suggesting runaway 
evaporation with increasing collision rate is possible. 

Time evolution of the phase space density <p(t) can be 
derived based on standard evaporation theory [16]. As- 
suming energetic atoms can leave the sample in all di- 
rections, we have (j)(t) = <ft(Q)(l + X^DOtftTot) 2 / and 
K3D = (jl — 5) / (77 — 4) [It], [3 • Here To is the initial colli- 
sion rate and X 3D = (ry-4)e-^/V2 [HI EH is the fraction 
of collisions producing an evaporated atoms. Here we see 
that a negative (3 < leads to a faster-than-exponential 
growth of the phase space density, which eventually di- 
verges at time t = (— Xsd^PTo) -1 . We compare the mod- 
els and our experiment result in Fig. 4. To reach the same 
final phase space density, the trap-tilting scheme would 
require a much shorter evaporation time than the weak- 
ening scheme. For comparison, a potential with fixed 
trap frequency (y = 0), e.g., radio- frequency based evap- 
oration in magnetic traps, permits an even stronger run- 
away effect, see Fig. 4. 

Remarkably, the performance of our evaporation is 
consistent with the 3D evaporation model. The con- 
sistency of our evaporation speed with the 3D model is 
somewhat surprising. In a strongly tilted trap where hot 
atoms can only escape the trap in the tilted direction, it is 
generally expected that the evaporation will exhibit per- 
formance consistent with one dimensional evaporation. 
In ID evaporation, energetic atoms are selected only by 
their velocity in the tilted direction. Integration over the 
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FIG. 5: (color online) Temperature and particle number de- 
pendence. Based on the 4 s evaporation data (solid circle), the 
temperature shows a polynomial dependence on particle num- 
ber T oc N a , with a = 1.46(2). For comparison, predictions 
from 3D (red) and ID (blue) evaporation models are shown 
using our experiment initial condition and r\ = 6.5 ± 0.3. 

Boltzmann distribution on only this degree of freedom 
then results in a reduction of the evaporation rate by a 
factor of 477 [HI and thus Aid = A3D/477. Performance of 
ID evaporation for all possible 77 is shown in the shaded 
area in Fig. 4. Our experiment result apparently permits 
evaporation speed much faster than the ID prediction. 

We suspect 3D-like evaporation in a tilted trap results 
from the inseparability of the potential and the existence 
of a saddle point located at the rim of the potential bar- 
rier. Both features, together with the similar trap fre- 
quencies in all directions, can lead to stochastic single 



particle motion [l9j]. When atoms with sufficiently high 
energy are created by collisions, stochastic motion allow 
them to efficiently find escape trajectories. If the ener- 
getic atoms have a high probability to escape, regardless 
of their initial direction of motion, evaporation is effec- 
tively three dimensional [l9j]. In realistic models, stochas- 
tization may also be induced by the intensity irregularity 
of the trapping laser beams. 

To further investigate the "dimension of evaporation" 
in a tilted trap, we come back to 77 + ft, which parame- 
terizes the energy removal by evaporating a single atom, 
or 77 + n = —(kBT)~ 1 dE/dN. For 3D evaporation, we 
expect k>sd = (?7 — 5)/(?7 — 4), which is ksd = 0.6(1) for 
our parameter 77 = 6.5(3); for ID evaporation, energy 
selectivity applies to the axial, but not the transverse 
motion, which has a mean energy of 2 ksT per particle. 
Hence, we expect a higher energy removal per particle 
with kid = &3D + 2 = 2.6(1) for our parameter [20[. Ex- 
perimentally, we can test these predictions by evaluating 
a = dlogT/dlogN, which has a simple dependence on k 
as shown in Eq. 4. We show in Fig. 5 that our 4 s evap- 
oration data is excellently fit to the polynomial function 
with a = 1.46(2). Using Eq. 4, we derive k = 0.6(3), 
which is consistent with the 3D value and confirms the 
3D nature of the trap-tilt based evaporation. 
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